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ABSTRACT 

We investigate the properties of the first galaxies at z > 10 with highly resolved 
numerical simulations, starting from cosmological initial conditions and taking into 
account all relevant primordial chemistry and cooling. A first galaxy is characterized 
by the onset of atomic hydrogen cooling, once the virial temperature exceeds ~ 10 4 K, 
and its ability to retain photoheated gas. We follow the complex accretion and star 
formation history of a ~ 5 x 10 7 M system by means of a detailed merger tree 
and derive an upper limit on the number of Population III (Pop III) stars formed 
prior to its assembly. We investigate the thermal and chemical evolution of infalling 
gas and find that partial ionization at temperatures > 10 4 K catalyses the formation 
of H 2 and hydrogen deuteride, allowing the gas to cool to the temperature of the 
cosmic microwave background. Depending on the strength of radiative and chemical 
feedback, primordial star formation might be dominated by intermediate-mass Pop III 
stars formed during the assembly of the first galaxies. Accretion on to the nascent 
galaxy begins with hot accretion, where gas is accreted directly from the intergalactic 
medium and shock-heated to the virial temperature, but is quickly accompanied by a 
phase of cold accretion, where the gas cools in filaments before flowing into the parent 
halo with high velocities. The latter drives supersonic turbulence at the centre of the 
galaxy and could lead to very efficient chemical mixing. The onset of turbulence in 
the first galaxies thus likely marks the transition to Pop II star formation. 
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INTRODUCTION 



Understanding the formation of the first stars and galaxies 
at the end of the cosmic dark ages is one of the most im- 
portant challenges in modern cosmology llBarkana fc Lo"eb 



20011 : iBromm fc Larsonll2004 ICiardi fc Ferrarall2005l ; lGlover 



nantly very massive, with M, > 100 Mp (lOmukai fc Palial 



2003l:lBromm fc Loeb|[20oi:IO'Shea fc Normanl2007l: but see 



Tan fc McKed 120041: iRipamontil 120071 ; iMcKee fc Tanl [20071 : 
McGreer fc Brvanll2008l ').~ 



20051 ). In the standard A cold dark matter (ACDM) cos- 
mology, the first stars, termed Population III (Pop III), 
are predicted to have formed in dark matter (DM) 
' minihaloes' with v ir ial mass ~ 10 6 Mq at z ~ 20 
(IBromm etTal] Il999l, |2002|; iNakamura fc Umemural l200ll; 
Abel et alj|200l |2002|; lYoshida et al.ll200l l200d : iGao et ail 
20071 : lO'Shea fc Normanll2007l ). Based on molecular hydro- 
gen cooling, the gas can cool to ~ 200 K and become Jeans- 
unstable once the central clump attains ~ 10 3 Mg. Such 
high temperatures lead to efficient accretion on to the pro- 
tostellar core and imply that the first stars were predomi- 
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Due to their primordial composition, massive Pop III 
stars have smaller radii than their present-day counterparts 
and their surface temperatures can exceed ~ 10 s K, 
resulting in h i gh ph ot on yields dTumlinson fc Shulll l2000l: 



Bromm et all 120011 ; ISchaererl l2002i : iTumlinson et al l 
20051 Consequently, they exert strong feedback on the 



intergalactic medium (IGM) via radiation in the Lyman- 
Werner (L W) bands, which rea di ly destroys mo l ecular 
hydrogen (I Glover fc Brand! l200ll: iMachacek et all 12001 



2006 



Ricotti et al.ll200ll: lOmukai fc Yoshiill2003l; iGreif fc Bromml 



MesinEcr et al 



Johnson et al 



2006 



■O'Shea fc Normanl 120081 : 
2007bl : IWise fc Abell l2007dl ). aiid ion- 
izing radiation giving rise to the first H 11 regions 
llKitavama et all 12004: IWhalen et all 120041: lAlvarez et all 



120061 : ISusa fc Umemural 120061 : lYoshida et ail 12007 
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Ahn fc Shapircl2007l : lAbel et al 1l2007l : I Johnson et al . 2007a; 



Whalen et al. 20071 ). Furthermore, a yet unknown fraction 



of massive, metal-free stars is expe cted to end in extremely 
violent supernova (SN) explosions l|Heger fc Wooslevl 1200a 



of dynamics and chemical enrichment 


Bromm et al. 2003; 


Mackev et al. 2003; 


Norman et al.l 20041; ISalvaterra et al.l 


20041; Yoshida et al. 


20041; Kitavama & Yoshidal 20051; 


Machida et all 12005; 


Greif et al.ll2007l; IWhalen et alj|2008h. 



Since primordial star formation initially occurs in mini- 
haloes that constitute the progenitors of the first galaxies, 
stellar feedback is expected to play an important role 
during their assembly. 

When and where did the first galaxies form? Ac- 
cording to the bottom-up character of structure forma- 
tion, as described by standard Press-Schechter theory 
l|Press fc Schechteij[l974T ). the first > 5 x 10 7 M Q haloes as- 
sembled at z > 10 via hierarchical merging. This mass scale 
is set by the onset of atomic hydrogen cooling once the virial 
temperature exceeds ~ 10 4 K; these objects have therefore 
frequently been termed 'atomic cooling haloes'. Henceforth, 
we will synonymously use the term 'first galaxy', as these 
haloes can retain photoheated gas and therefore might, 
for the first time, maintain self-regula ted star formation 



19961: iMac Low & Ferraral 1 19991; 


Madau et al. 200 ll; 


Mori et al. 


20021; Oh & Haimanl 


20021; Ricotti et all 


2002313; 1 


Scannapieco et al. 2002; Thacker et al. 2002; 


Wada & Venkatesanl 20031; iDiikstra et al.ll2004l; iRead et all 


20061; Ricotti et al. 2008). 



What types of stars are expected to form during the 
assembly of the first galaxies? Depending on the primor- 
dial initial mass function (IMF) and the strength of ra- 
diative and SN-driven feedback, the IGM could be heav- 
ily enriched with metals prior to the onset of second- 
generation star formation (|Bromm et al.l 120031 ; iGreif et al.l 
l2007l ; IWise fc Abelll2007cl ). The existence and precise value 
of a critical metallicity for the transition to Pop II 
are still vigorously debated, and has originally been dis- 
cussed in the context of fine-structure versus dust cooling 



1 Bromm fc Loebll2003bl; ISchneider et al 



_ -..20031; iFrebel et all 

l2007l ; |jappsen et al.ll2007l ; ISmith fc Sigurdssonll2007l ) . In the 
former case, critical metallicities of Z cr it — 10 Zq have 



been found (|Bromm et al.l l200ll ; ISantoro fc Shulll [2006). 
while in the latter case uncertainties in the dust compo- 
sition and gas-phase depletion lead t o a range of possi- 
ble values, 10"^ < Zc rit < 10~ 5 Z p JOmukai et all |2005| ; 
ISchneider et~aill2006l ; iTsuribe fc Omukaill2006l ). Somewhat 
transcending this debate, recent simulations have shown 
that a single SN explosion might en rich the local IGM 
to well above any critical metallicit y (|Bromm et alj 120031 ; 
IGreif et al.ll2007l ; IWise fc Abelll2007d) . forcing an early tran- 
sition in star formation mode (|Clark et al.| [2008). However, 
if radiative feedback was sufficiently strong, the bulk of 
primordial star formation might ha ve occurred in systems 
that cool via atomic hydrogen lines llGreif fc Bromml 120061 ; 
lO'Shea fc Normanl [20081 : 1 Johnson et al.ll2007bl ). irilight of 
these uncertainties, we investigate the formation of a first 
galaxy in the limiting case of no feedback, allowing us to fo- 
cus on the chemistry, cooling and development of turbulence 
during the assembly of the atomic cooling halo. 

Theoretical investigations have recently pointed to- 



wards the existenc e of two physically distinct populations 
of me t al-free stars llUehara fc Inut suka 2000j; iMackev et al.l 
120031 ; Ijohnson fc Bromml 120061 )! Gas cooling primarily 
via molecular hydrogen leads to the formation of > 
100 Mq stars, termed Pop III.l, while in regions of 
previous ionization hydrogen deuteride (HD) likely en- 
ables the formation of > 10 Mq stars, termed Pop III. 2 
jTan fc McKed I2OO8I ). This latter mode had previously 
been termed Pop II.5 ilMackev et al.|[2003l ; IGreif fc Bromml 
120061 ; Ijohnson fc Bromml 120061 ). Scenarios providing a sub- 
stantial degree of ionizat ion include relic H n regions 
l|Nagakura fc Om ukai 20 051), dense shel l s produced by ener- 
getic SN explosions (iBromm et al.ll20q3l;lMackev et al]|2003l ; 
ISalvaterra et al.l |2004| ; iMachida et al.ll2005l ). and structure 
formation shocks in the virialization of th e first galaxies 
jOh fc Haimanl 120021 ; [Greif fc Brom m 2006). In the former 
case, numerical simulations have confirmed that the gas can 
cool to the temperature of the cosmic microwave background 
(CMB), reducing the Jeans- mass by almost a n order of mag- 
nitude compared to the tru ly primordial case l|johnson et al.l 
l2007al ; lYoshida et al.|[2007l) , while we here set out to investi- 
gate the effects of partial ionization and molecule formation 
during the assembly of the first galaxies. One of the key 
questions is therefore whether the gas can cool to the CMB 
limit and thus enable the formation of Pop III. 2 stars. 

Another important aspect concerning the baryonic col- 
lapse of the first galaxies is the development of turbulence. 
While turbulence does not seem to play a role in mini- 
haloes, at least on s cales comparable to their virial radius 
|Yoshida et al.| [2006). turbulent motions could become im- 
portant in more massive haloes, where accretion of cold gas 
from filaments in combinatio n with a softened equation of 
state drives strong shocks l|Keres et alj 120051 ; IWise fc Abell 
l2007bl ; ISancisi et al.l l2008h . possibly leading to vigorous 
fragmentation and the formation of the first star clusters 
l|Clark et al.ll2008l ). This would mark the first step towards 
conditions similar to present-day star formation, where su- 
personic turbulent velocity fields determine the fragmenta- 
tion properties of the gas. In our present study, we inves- 
tigate the role of turbulence by analysing the velocity field 
and energy content of the galaxy during virialization. We 
then discuss the likely fragmentation properties of the gas 
and the consequences for second-generation star formation. 

The final ingredient of early galaxy formation is the 
co-evolution of massive black holes (MBHs) and the sur- 
rounding stellar system, leading in extreme cases to the 
formation of the highest redshift quasars observed at 
z > 6. The most sophisticated simulations of supermas- 
sive black hole (SMBH) growth to date have been con- 
ducted by iLi et all (|2007T ). but their resolution was insuf- 
ficient to trace the origin of SMBHs below ~ 10 5 M . 
Such SMBHs could form via direct collapse of isother- 
mal gas in atomic cooling haloes in t he presence of 
a stro n g photo-dissociatin g back g round (jBromm fc Loebl 
2003al; iKoushiappas et al.l 1200-3; iBegelman et al.l 120*061 



Lodato fc NataraiarJl2006tlSpaans fc SiUdl2006t IWise et all 
20071 ). or via merging and accreti on of gas on t o the 



Schneider et al. 20021; 


Islam et al.l 2003; Volonteri et al.l 


2003; Volonteri & Rees 


2005; 


Johnson & Bromml 120071: 


Pelupessv et al.ll2007l; Wise et al.l 


2007). Here, we investigate 



the fate of BHs seeded by the remnants of Pop III stars and 
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their prospect of growing into SMBHs at the centres of the 
first galaxies. Finally, we provide an estimate for the ampli- 
tude of ionizing and molecule-dissociating radiation emitted 
by the accretion disc around the central BH. 

The structure of our work is as follows. In Section 2, we 
discuss the setup of the cosmological simulations, the im- 
plementation of multiple levels of refinement and our sink 
particle algorithm. We then describe the hierarchical assem- 
bly of the galaxy (Section 3), its cooling and star formation 
properties (Section 4) and the development of turbulence 
(Section 5). In Section 6, we discuss the growth of the BH 
at the centre of the galaxy, and in Section 7 we summarize 
our results and assess their implications. For consistency, all 
distances are physical (proper), unless noted otherwise. 



2 NUMERICAL METHODOLOGY 

To investigate the formation of the first galaxies with ad- 
equate accuracy, our simulations must resolve substructure 
on small scales as well as tidal torques and global gravi- 
tational collapse on much larger scales. We capture all the 
relevant dynamics by performing a preliminary simulation 
with a coarse base resolution, but refine around the point 
of highest fluctuation power. This ensures that the region 
containing the mass of the galaxy is well resolved and that 
a given number of particles is distributed efficiently. 

2.1 Initial setup 

The simulations are carried out in a cosmological box of 
linear size ~ 700 kpc (comoving), and are initialized at 
z = 99 according to a concordance A cold dark matter 
(ACDM) cosmology with matter density O m = 1 — Aa = 
0.3, baryon density f2b = 0.04, Hubble parameter h = 
H / (100 km s _1 Mpc" 1 ) = 0.7, spectral index Tig = 1.0, 
and normalization as — 0.9 ( Sp ergel et alj [20031 ) . Density 
and velocity perturbations are imprinted at recombination 
with a Gaussian distribution, and we apply the Zeldovich ap- 
proximation to propagate the fluctuations to z = 99, when 
the simulation is started. To capture the chemical evolu- 
tion of the gas, we follow the abundances of H, H + , H~, 
H2, H^~, He, He + , He ++ , and e~, as well as the five deu- 
terium species D, D + , D~, HD and HD + . We include all rel- 
evant cooling mechanisms, i.e. H and He atomic line cooling, 
bremsstrahlung, inverse Compton scatte ring, and collisional 
excit ation cooling via H2 and HD (see iJohnson fe Bromml 

l2ooeh . 



of 512 child particles. The particle mass in the region con- 
taining the comoving volume of the galaxy with an extent 
of ~ 200 kpc (comoving) is ~ 100 M Q in DM and ~ 10 M Q 
in gas. This allows a baryonic resolution of ~ 10 3 M@, such 
that the 'loitering state' in minihaloe s at T ~ 200 K an d 
nu — 10 4 cm -3 is marginally resolved (|Bromm et alj|2~002h . 

To follow the growth of BHs seeded by the collapse 
of Pop III. 1 stars in minihaloes, we apply a slightly mod- 
ified version of the sink particle algorithm introduced in 
IJohnson fc Bromml (|2007T ). Specifically, we create sink par- 
ticles once the density exceeds riH = 10 4 cm -3 and immedi- 
ately accrete all particles within the resolution limit of the 
highest-density particle, resulting in an initial sink particle 
mass of ~ 2 x 10 3 Mq . Further accretion is governed by the 
following criteria: particles must be bound (i?kin < |-E po t|), 
have a negative divergence in velocity and fall below the 
Bondi radius: 

fimnGM B H /lN 

where /i is the mean molecular weight, Mbh the mass of the 
black hole, and T the temperature of the gas at the location 
of the sink particle, determined by a mass-weighted aver- 
age over all accreted particles. Typical values for the initial 
Bondi radius are re — 5 pc. Here, we do not model radiative 
feedback that would result from accretion, and we therefore 
overestimate the Bondi radius and consequently the growth 
of the BH. We further assume that BH mergers occur in- 
stantaneously once their separation falls below the resolu- 
tion limit, provided that their relative velocity is smaller 
than the escape velocity. Our idealized treatment provides a 
robust upper limit on the growth rate of BHs during the as- 
sembly of the first galaxies, but a more realistic calculation 
including the effects of feedback will be presented in future 
work. 



3 FIRST GALAXY ASSEMBLY 

Among the key questions concerning the formation of the 
first galaxies are the degree of complexity associated with 
the halo assembly process, the role of previous star forma- 
tion in minihaloes, the chemical and thermal evolution of 
infalling gas, and the development of turbulence. To answer 
these questions, we first discuss the gravitational evolution 
of the DM with a merger tree that reconstructs the mass 
accretion history of the resulting atomic cooling halo, and 
allows us to determine the maximum amount of previous 
star formation in its progenitor haloes. 



2.2 Refinement and sink particle formation 

In a preliminary run with 64 3 particles per species (DM and 
gas) , we locate the formation site of the first ~ 5 x 10 7 Mq 
halo. This object is just massive enough to activate atomic 
hydrogen cooling and fulfil our prescription for a galaxy. 
We subsequently carry out a standard hierarchical zoom-in 
procedure to achieve high mass resolutio n inside the region 
destined to collapse into the galaxy fe.g. [Navarro fe White! 
1 19941 : iTormen et ai1ll997l ; iGao et al.ll2005l ). We apply three 
consecutive levels of refinement centred on this location, 
such that a single parent particle is replaced by a maximum 



3.1 Atomic cooling criterion 

Possibly the most important characteristic distinguishing 
the first galaxies from their lower-mass minihalo predeces- 
sors is their ability to cool via atomic hydrogen lines, which 
softens the equation of state at the virial radius and allows a 
fraction of the potential energy to be converted into kinetic 
energy. At high redshift, the virial temperature of a system 
with virial mass M V [ T can be expressed as 

— 4K (OTf(T)< m 
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such that a 5 x 10 7 Ma halo at 



10 i s just massive 



enou gh to fulfil the atomic cooling criterion (|Oh fc Haimanl 
2002). Related to this, they have the ab ility to retain 
gas photoheated by hydrogen ionization (|Diikstra et al.l 
120041 '). likely allowing self-regulated star formation for the 
first time. Furthermore, it is often argued that star for- 
mation in atomic cooling haloes provided the bulk of 
the photons for r eionization due to e f ficient shielding 
from L W radiation |Bromm fc Loeb 2003a; Greif fc Bromml 



Press-Schechter 
J Simulation 



20061; iHaiman fc Brvanl 120061 ; IChoudhurv et all l200sf : 
Wise fc Abelll2007al ). The y may have been t he key drivers 
of early IGM enrichment (|Madau et al.ll2001T l. and the host 
systems for the formation of the first low-mass stars that 
can be pro bed in the present-day Milky Way via st ellar ar- 
chaeology (|Frebel et al.ll2007l ; Kaxlsson et al.ll2007l ). Their 
formation thus marks an important milestone in cosmic his- 
tory, and in light of upcoming observations it is particularly 
important to understand the properties of these systems. 



3.2 Cosmological abundance 

How common are atomic cooling haloes at the redshifts we 
consider? A frequently used tool to estima te their abundance 
is th e Press-Schechter mass function (|Press fc Schechterl 
1974), which provides an analytic expression for the num- 
ber of DM haloes per mass and comoving volume. Even 
though the Press-Schechter mass function is inaccurate at 
low redshifts compared to other a nalytic estimates such a s 
the Sheth-Tormen mass function (|Sheth fc Tormer] 2002), 
it provides a good fit to n umerical simulations at the early - 
times we are interested in (Ijang-Condell fc Hernquist|[200ll ; 
iHeitmann et al.ll2006l ; iReed et al.ll2007T l. In Fig. 1. we show 
the Press-Schechter mass functions for z ~ 30 and 10 (solid 
and dotted line, respectively), while the symbols indicate 
the halo distribution extracte d from the simulation b y the 
group-finding algorithm HOP l|Eisenstein fc Huuli"998T ). The 
simulation and the analytic estimate generally agree, al- 
though there is a large scatter due to the finite box size. 
Better agreement could be expected if one accounted for 
the fluctuation power on scales gr eater than the computa- 
tional box (see lYoshida et aLll2003l ). At redshift z ~ 10, one 
expects to find roughly 10 atomic cooling haloes per cubic 
Mpc (comoving), i.e. of the order of unity in our computa- 
tional box of length ~ 700 kpc (comoving). Indeed, in our 
simulation a single 5 x 10 7 Mq halo forms by z ~ 10. 



3.3 Assembly of atomic cooling halo 

In the ACDM paradigm, structure formation proceeds hi- 
erarchically, with small objects collapsing first and sub- 
sequently growing via merging. This behaviour eventually 
leads to the formation of haloes massive enough to fulfil 
the atomic cooling criterion. In Fig. 2, we show the DM 
overdensity, gas density and temperature averaged along the 
line of sight at three different output times. The brightest 
regions in the DM distribution mark haloes in virial equilib- 
rium, according to the commonly used criterion p/p ~ 178, 
where p and p denote the local and background density. 
White crosses denote the formation sites of Pop III. 1 stars 
in minihaloes. The first star-forming minihalo at the cen- 
tre of the box assembles at z ~ 23 and subsequently grows 




Figure 1. Comparison of the Press-Schechter mass function 
(solid and dotted line) and the simulation results (crosses and 
triangles) at z ~ 30 (lower set) and z ~ 10 (upper set). The sim- 
ulation generally agrees with the analytic prediction, although 
there is a large scatter due to the finite box size. The expected 
number of atomic cooling haloes per cubic Mpc (comoving) at 
z ~ 10 is of the order of 10. 



into the galaxy delineated by the insets in the right panels 
of Fig. 2 and further enlarged in Figs 8, 9, 10 and 12. Al- 
though this structure is not yet fully virialized and exhibits 
a number of sub-components, it has a common potential well 
and attracts gas from the IGM towards its centre of mass, 
where it is accreted by the central BH once it falls below the 
Bondi radius. The virial temperature of the first minihalo in- 
creases according to equation (2) until it reaches ~ 10 4 K 
and atomic cooling sets in, at which point the equation of 
state softens and a fraction of the potential energy is con- 
verted into kinetic energy. Star formation takes place only 
in the most massive minihaloes, with 10 Pop III. 1 star for- 
mation sites residing in the volume that is destined to form 
the galaxy. This has important consequences for the role of 
stellar feedback, and will be further discussed in Section 4. 
The morphology of the galaxy can best be seen in Fig. 3, 
where we show a three-dimensional rendering of the central 
150 kpc (comoving), i.e. the same field of view as in Fig. 2. 
The temperature is colour-coded such that the hottest re- 
gions with T ~ 10 4 K are displayed in bright red. Here, the 
true spacial structure of the galaxy becomes more apparent, 
showing that its environment is organized into prominent 
filaments with a high amount of substructure. In some in- 
stances, star-forming minihaloes have aligned along these 
filaments and will soon merge with the galaxy. 



3.4 Merger tree 

The hierarchical assembly of the galaxy can be best de- 
scribed by means of a merger tree that depicts the evolution 
of all progenitor haloes. We construct such a merger tree by 
tagging all DM particles that reside in the parent atomic 
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Figure 2. The DM overdensity, hydrogen number density and temperature averaged along the line of sight within the central ~ 150 kpc 
(comoving) at three different output times, from z ~ 23, when the first star-forming minihalo at the centre of the box collapses, to 
z ~ 10, when the first galaxy forms. White crosses denote Pop III. 1 star formation sites in minihaloes, and the insets approximately 
delineate the boundary of the galaxy, further enlarged in Figs 8, 9, 10 and 12. Top row: The hierarchical merging of DM haloes leads 
to the collapse of increasingly massive structures, with the least massive progenitors forming at the resolution limit of ~ 10 4 Mq and 
ultimately merging into the first galaxy with ~ 5 X 10 7 Mq. The brightest regions mark haloes in virial equilibrium according to the 
commonly used criterion p/p > 178. Although the resulting galaxy is not yet fully virialized and is still broken up into a number of 
sub-components, it shares a common potential well and the infalling gas is attracted towards its centre of mass. Middle row: The gas 
generally follows the potential set by the DM, but pressure forces prevent collapse in haloes below ~ 2 X 10 4 Mq (cosmological Jeans 
criterion). Moreover, star formation only occurs in haloes with virial masses above ~ 10 5 Mq, as densities must become high enough for 
molecule formation and cooling. Bottom row: The virial temperature of the first star-forming minihalo gradually increases from ~ 10 3 K 
to ~ 10 4 K, at which point atomic cooling sets in. 



cooling halo and track their location backwards in time. If 
they are part of a group at a previous timestep, they are 
considered to reside in a halo with mass equal to the sum 
of their individual masses. We repeat this process until all 
tagged particles are no longer part of a group or the mass 
falls below the resolution limit. Fig. 4 shows the resulting 
absolute and differential mass growth of the galaxy. Accre- 



tion is fuelled by minor as well as major mergers, with the 
latter showing the tendency to double the mass of the halo. 
In the course of ~ 400 Myr, the accretion rate increases 
from ~ 5 x 10" 3 M yr" 1 to ~ 0.5 M Q yr -1 , but varies 
significantly in between due to the highly complex nature of 
bottom-up structure formation. 

To illustrate the substantial degree of complexity in- 
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Figure 3. A three-dimensional rendering of the central ~ 150 kpc (comoving), showing the same field of view as in Fig. 2. The temperature 
is colour-coded such that the hottest regions with T ~ 10 4 K are displayed in bright red. Here, the true spacial structure of the galaxy 
becomes more clear, showing that its environment is organized into prominent filaments with a high amount of substructure. In some 
instances, star-forming minihaloes have aligned along these filaments and will soon merge with the galaxy. 



volved, Fig. 5 shows the individual paths of all progeni- 
tor haloes down to the DM resolution limit of ~ 10 4 Mq. 
The initial widening of the tree indicates that an increas- 
ing number of minihaloes collapse, while at z ~ 20 merg- 
ing becomes dominant and the degree of complexity de- 
creases again. The timescale for the widening of the tree 
is ~ 150 Myr, while the completion of the merging process 
requires another ~ 250 Myr. The total number of haloes 
above ~ 10 4 Mq that merge to form the galaxy is ~ 300. 
A more sophisticated analysis of the merging process is pre- 
sented in Fig. 6, where the individual paths and masses of 
all progenitor haloes are shown. Each line represents an in- 
dividual halo, while the colour denotes its mass. The target 
halo seeding the first galaxy is indicated by the rightmost 



path. Sites of Pop III. 1 star formation are denoted by star 
symbols, and the oval denotes the formation of the atomic 
cooling halo. The history of this most basic building block 
of galaxy formation is highly complex, further complicated 
by the formation of 10 Pop III. 1 stars prior to its assembly. 
However, this is an upper limit on previous star formation 
activity as we do not include radiative and SN-driven feed- 
back, which would likely reduce the net star formation rate. 

The presence of DM fluctuations on mass scales be- 
low our resolution limit might imply that Pop III star for- 
mation takes place in haloes with viral mass well below 
~ 2 x 10 4 Mq , but pressure forces prevent gas from settling 
into these shallow potential wells (cosmological Jeans crite- 
rion). Moreover, gas may not be able to collapse beyond the 
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t H [Myr~ 
128 177 




Figure 4. The virial mass (top panel) and accretion rate (bottom 
panel) of the galaxy as a function of redshift. The growth of the 
underlying DM halo is fuelled by minor as well as major mergers, 
with the latter showing the tendency to double the mass of the 
target halo. At z ~ 10, the atomic cooling criterion is fulfilled and 
a galaxy is born. 



point of virialization in < 10 5 Mq haloes, as temperatures 
and densities do not become high enough for efficient H2 
formation. Dynamical heating by mergers counteracts cool- 
ing and thus only a fraction of a ll minihaloes will be able 
to form stars l|Yoshida et al.l l2003). To ascertain the impor- 
tance of this effect, we determine the masses of all minihaloes 
experiencing star formation. As described in Section 2, the 
formation of a Pop HIT star is denoted by the creation of 
a sink particle once the hydrogen number density exceeds 
~ 10 4 cm~ 3 . We implicitly assume that such a parcel of 
gas does not experience further subfragmentation. With this 
prescription, we find the following virial masses for all star- 
forming minihaloes shown in Fig. 6, from top to bottom: 
M vir ~ [5.8, 1.6, 7.5, 3.5, 4.3, 1.4, 3.2, 9.3, 1.4, 11.8] x 10 5 Mq. 
As expected, their masses are in the range ~ 10 s - 10 6 Mq, 
emphasizing the influence of dynamical heating on haloes 
below ~ 10 5 Mq . Interestingly, the fact that only a fraction 
of all minihaloes forms stars ensures a constant inflow of 
cold gas into existing haloes, which is crucial for the growth 
of the BHs at their centres. 



4 COOLING AND STAR FORMATION 

A crucial issue concerning the formation of the first galaxies 
is the chemical and thermal evolution of accreted gas, which 
ultimately determines the mode of star formation. We here 
briefly discuss radiative and SN-driven feedback exerted by 
the very first stars, followed by a discussion of the chemistry 
and cooling properties of an atomic cooling halo and the 
implications for second-generation star formation. 



U z 



"*■ > 

1 11 1. > 




79 



Figure 5. The merger tree of the galaxy, illustrating its complex- 
ity as a function of time. New branches indicate the formation of 
DM haloes at the resolution limit of ~ 10 4 Mq. The widening of 
the tree increases as more and more haloes collapse, until merg- 
ing dominates and the degree of complexity decreases again. The 
timescale for the former is ~ 150 Myr, while the completion of 
the merging process requires another ~ 250 Myr. A total of ~ 300 
haloes above ~ 10 4 Mq merge to form the galaxy. 



4.1 Population III.l 

As shown in Section 3, star formation ensues in mini- 
haloes before the larger potential wells of the first galax- 
ies assemble. This implies that radiative and SN-driven 
feedback influences star formation in other minihaloes as 
well as second-generation star formation in the resulting 
atomic cooling halo. However, recent numerical simulations 
have shown that local radiative feedback via photoheat- 
ing and LW ra diation may not be as important as previ- 
ously thought (lAhn fc Shapiro! 120071 ; I Johnson et"aH l2007al ; 
IWhalen et al.ll2007 ). and that a global LW background may 
only reduce the number of Pop HIT stars by ~ 50 per cent 
jGreif fc Brommll2006l ; Ijohnson et al.ll2007bi ). An unknown 
fraction of these stars end their lives as energetic SNe and 
enrich the surrounding IGM to well above the critical metal- 
licity (|Bromm et al.l 120031 ; iGreif et all 120071 ; IWise fc Abel 
l2007cl ). while ot hers collapse directly to BHs and do not ex - 
pel any metals (|Heger fc Wooslevl |2002| ; iHeger et al.ll2003l ). 
Since the t imescale for the recollapse of enriched gas is 
> 100 Myr l|Greif et al.ll2007l) . and mixing is inefficient with 
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Figure 6. The full merger tree of the galaxy assembling at z ~ 10. Each line represents an individual progenitor halo and is colour-coded 
according to its mass. The target halo seeding the galaxy is represented by the rightmost path, which ultimately attains ~ 5 X 10 7 Mq 
and fulfils the atomic cooling criterion (denoted by the red oval). Star symbols denote the formation of Pop III. 1 stars in minihaloes, 
showing that in our specific realization 10 Pop III. 1 stars form prior to the assembly of the galaxy. Only a fraction of all minihaloes form 
stars, as dynamical heating via mergers partially offsets cooling. Depending on the detailed merger history, this ensures that star-forming 
minihaloes are supplied with cold gas, which is crucial for the growth of the BHs at their centres. 
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Figure 7. The phase-space distribution of gas inside the first star-forming minihalo (left-hand panel) and the atomic cooling halo 
(right-hand panel). We show the temperature, electron fraction, HD fraction and H2 fraction as a function of hydrogen number density, 
clockwise from top left to bottom left. Left-hand panel: In the minihalo case, adiabatic collapse drives the temperature to > 10 3 K and 
the density to njj > 1 cm -3 , where molecule formation sets in allowing the gas to cool to ~ 200 K. At this point, the central clump 
becomes Jeans-unstable and ultimately forms a Pop III. 1 star. Right-hand panel: In the first galaxy, a second cooling channel has emerged 
due to an elevated electron fraction at the virial shock, which in turn enhances molecule formation and allows the gas to cool to the 
temperature of the CMB. The dashed red lines and arrows approximately delineate the resulting Pop III. 1 and Pop III. 2 channels, while 
the solid green lines denote the path of a representative fluid element that follows the Pop III. 2 channel. 



respe ct to pre-established overdensities (|Cen fc Riquelmel 
120081 ). subsequent star formation in minihaloes prior to the 
assembly of the atomic cooling halo likely remains metal- 
free. It is much more difficult, however, to predict the char- 
acter of star formation inside the first galaxies. In the fol- 
lowing, we will first examine the consequences of pristine 
gas collapsing in the atomic cooling halo, and subsequently 
briefly address the corresponding case of pre-enriched gas. 



4.2 Population III. 2 

The possible existence of a distinct population of metal- 
free stars in regions of previous ionizati o n has attracted 



increasing attention (IMackev et al.l 12003; iGreif & Bromm 


2006 


;l Johnson & Brommll2006l;lYoshidall2006l: Tan & McKee 


2008 


; Yoshida et al.l 120071; but see iRipamonti 2007; 


McGreer & Brvan 2008). According to theory, an elevated 



electron fraction catalyses the formation of H2 and HD well 
above the level found in minihaloes and enables the gas 
to cool to the temperature of the CMB. This reduces the 
Bonnor-Ebert mass by almost an order of magnitude and 
likely leads to the formatio n of Pop III. 2 stars with > 10 Mq 
jjohnson fc Bromml 120061 ). Numerical simulations of star 
formatio n in relic H 11 regions have largely confirm ed this 
picture (jjohnson et all l2007al ; lYoshida et all l2007h . while 
its relevance during the virialization of the first galaxies 



has not yet been esta blished (but see IGreif fc Brommll2006l ; 
iJohnson et"aill2007bl ). 

The chemistry of gas contracting in an atomic cooling 
halo is fundamentally different from that in minihaloes. The 
latter maintain a primordial electron fraction of ~ 3 x 10~ 4 
and form a limited amount of molecules, while the virial tem- 
perature in an atomic cooling halo exceeds ~ 10 4 K and the 
elevated electron fraction facilitates the formation of high 
H2 and HD abundances. This allows the gas to cool to the 
temperature of the CMB instead of the canonical ~ 200 K 
found in minihaloes. In the left panel of Fig. 7, we show the 
properties of the gas in the first star-forming minihalo at 
the centre of the computational box. The primordial elec- 
tron fraction remains constant until densities become high 
enough for electron recombination. After adiabatic heating 
to > 10 3 K, molecule formation sets in and the gas cools to 
~ 200 K, at which point the central clump becomes Jeans- 
unstable and inevitably forms a Pop III. 1 star. In contrast, 
the right panel of Fig. 7 shows the density and temperature 
of the gas inside the atomic cooling halo. The conventional 
H2 cooling channel is still visible, but a second path from 
low to high density has emerged, enabled by an elevated 
electron fraction at the virial shock, which in turn enhances 
the formation of H2 and HD and allows the gas to cool to the 
temperature of the CMB. The dashed red lines and black ar- 
rows in Fig. 7 approximately delineate both channels, show- 
ing that the electron fraction in the Pop III. 2 case is elevated 
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by an order of magnitude to ~ 10 3 , and t he H2 fraction 
rises to ~ 2 x 10 -3 . As already estimated in I Johnson et all 
11200701 ) , the HD fraction grows to above ~ 10" 6 . 

To more clearly illustrate this point, we plot the path of 
a representative fluid element evolving along the Pop III. 2 
channel (solid green lines in Fig. 7). Such gas indeed cools 
to the CMB floor, potentially enabling the formation of 
Pop III. 2 stars. The mass fraction entering this channel 
is relatively low at the time considered here since the 
atomic cooling threshold has just been surpassed, but should 
quickly rise as freshly accreted material is shock-heated to 
~ 10 4 K. We cannot study any possible fragmentation, since 
the gas rapidly falls to within the Bondi radius of the cen- 
tral BH, and is thus accreted by the sink particle. It will 
be very interesting to investigate the fragmentation of the 
Pop III. 2 mode in future, higher-resolution simulations, in 
par ticular testing the predicted mass scale of > 10 Mq (see 
e.g. IClark et al.ll2008l ). 

What are the implications of this result? Parcels of 
gas that are accreted on to the galaxy through the virial 
shock can cool to the temperature of the CMB and pos- 
sibly become gravitationally unstable, resulting in the for- 
mation of Pop III. 2 stars. Including radiative feedback from 
previous star formation would only strengthen this conclu- 
sion, as the degree of ioniz ation would be increased even 
further (|Wise fc Abelll2007cD . As long as the gas collapsing 
into the first galaxies remains pristine, primordial star for- 
mation will therefore li kely be dominated by intermediate- 
mass (Pop III. 2) stars (|Greif fc Brommll2006r ). The crucial 
question, however, is: Can the gas inside the first galaxies 
remain metal-free? 



4.3 Population II 

In the previous sections, we have found that of order 10 
Pop III. 1 stars form prior to the assembly of the atomic 
cooling halo. In this case it appears unlikely that all of 
them will collapse into BHs without any metal-enrichment 
(| Johnson et al.l T2007bl ). Even a single SN from a massive 
Pop III star would already suffice to reach levels above 
the critical metallicity. at least on average (|Bromm et al.l 
120031 ; iGreif et all 120071 ; IWise fc Abel l2007d ). More gener- 
ally, at some stage in cosmic history, there must have been 
a transition from primordial, high-mass star formation to 
the 'normal' mode that dominates today. The discovery of 
extremely metal-poor st ars in the Galactic hal o with masses 
below one solar mass dChristlieb et al l |2002| ; | Frebel et al.l 
l2005l ; lBeers fc Christliebll2005f ) indicates that this transition 
occurs at abundances considerably smaller than the solar 
value. At the extreme end, these stars have iron abundances 
less than 10 -5 times the solar value, but show significant 
carbon and oxygen enhancements, which could be due to 
unusual abundance patter ns produced by enrichm ent from 
BH-forming Pop III SNe (|Umeda fc Nomotd [2001 ) . or due 
to mass transfer from a close binary companion, whose fre- 
quency is predicted t o increase with decreasing metallicity 
l|Lucatello et alj|2005h . 

Identifying the critical metallicity at which this transi- 
tion occurs is subject to ongoing research. One approach 
is to argue that low mass star formation becomes pos- 
sible only when atomic fine-structur e line cooling from 
carbon and oxygen becomes effective (jBromm et al. l200ll ; 



Bromm fc Loebl l2003bl ; ISantoro fc Shun1l2006l ; iFrebel et all 
20071 ). setting a value for Z c 



L cr jt at 



10" 



Another 



possibility, first proposed by lOmukai et all (|2005T ). is that 
low mass star formation is a result of dust-induced fragmen- 
tation occurring at high densities, iih — 10 13 cm" 3 , and 
thus at a very late stage in the protostellar collapse. In this 



model, 10 



< Zcrit < 10- 



where much of the un- 



certainty in the predicted value results from uncertainties 
in th e dust composition and the degree of gas-phase deple- 
tion dSchneider et alj|2002l. 120061). Recent numerical simu- 
lation s by iTsuribe fc Omukail (|20o3 ) as well as lClark et al.l 
(2008) provide support for this picture. However, the exist- 
ing data of metal-poor Galactic halo stars seems to be well 
accommodated by the C- and O-based fine-structure model 
(|Frebel et al.ll2007l ). 

In the present simulation, we do not follow the metal- 
licity evolution of the infalling gas. Thus, we can only specu- 
late about the properties of the resulting stellar population. 
It appears reasonable to assume that some of the accret- 
ing material is still pristine and free of metals, triggering 
the formation of lower-mass metal-free Pop III. 2 stars. Gas 
that flows in at even later times may already have experi- 
enced metal enrichment from previous Pop III SNe in nearby 
minihaloes. Because of the high level of turbulence within 
the virial radius at that time (see Section 5), the incoming 
new material is likely to efficiently mix with the pre-existing 
zero-metallicity gas and the era of Pop III star formation 
could come to an end. This transition possibly occurs at 
the same time as the onset of significant degrees of turbu- 
lence in the atomic cooling halo. We therefore speculate that 
some of the extremely metal-deficient stars in the halo of 
the Milky Way ma y have formed as early as redshift z ~ 10 
(|Clark et al.ll2008l ). 



5 TURBULENCE 

The development of turbulence in gas flowing into the cen- 
tral potential well of the halo strongly influences its fragmen- 
tation behaviour and consequently its ability to form stars. 
Detailed studies of the interstellar medium in the Milky 
Way, for instance, tell us that turbulence determines when 
and where star formation occurs and that it is the intricate 
interplay between gravity on the one hand, and turbulence, 
thermal pressure and magnetic fields on the o ther that sets 
the properties of young stars and star clusters ijLarsonl 120031 ; 
iMac Low fc Klessen|[200i ; iBallesteros-Paredes et al.l 120071 ) . 
In the context of our work, we investigate the velocity field 
and energy distribution that build up during the assembly of 
the galaxy. As opposed to cooling flows in low-mass haloes, 
the accretion flow on to the deep central potential well of 
the atomic cooling halo considered here becomes highly tur- 
bulent within the virial radius. 



5.1 The development of turbulence: hot versus 
cold accretion 

One of the most important consequences of atomic cooling 
is the softening of the equation of state below the virial 
radius, allowing a fraction o f the potential energy to be con- 
verted into kinetic energy (|Wise fc Abell l2007b! ) . This im- 
plies that perturbations in the gravitational potential can 
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generate turbulent motions on galactic scales, which are then 
transported to the centre of the galaxy (e.g. Fig. 8). In this 
context it is important to investigate the accretion of gas on 
to the galaxy in more detail. 

In principle, there are two distinct modes of accretion. 
Gas accreted directly from the IGM is heated to the virial 
temperature and comprises the sole channel of inflow un- 
til cooling in filaments becomes important. This mode is 
termed hot accretion, and dominates in low-mass haloes at 
high redshift. In the atomic cooling halo, the formation of 
the virial shock and the concomitant heating are visible in 
Fig. 9, where we show the hydrogen number density and 
temperature of the central ~ 40 kpc (comoving) around the 
BH at the centre of the galaxy. This case also reveals a sec- 
ond mode, termed cold accretion. It becomes important as 
soon as filaments are massive enough to enable molecule 
reformation, which allows the gas to cool and flow into the 
central regions of the nascent galaxy with high velocities. In 
Fig. 9, the cold gas accreted along filaments from the left- 
and right-hand is clearly distinguishable from the hot gas at 
the virial shock. These streams are also visible in Fig. 10, 
where we compare the radial with the tangential velocity 
component, and in Fig. 8, where we show the Mach number 
of infalling gas. Evidently, inflow velocities can be as high 
as 20 km s _1 , with Mach numbers of the order of 10. 

In Fig. 11, we compare the energy distribution and mass 
fraction of cold (< 500 K) versus hot (> 500 K) gas in radial 
shells for the first star-forming minihalo just before the for- 
mation of the sink particle, and the atomic cooling halo as- 
sembling at z ~ 10. The blue, green and red lines denote the 
azimuthally averaged ratio of radial, tangential and thermal 
to potential energy, respectively. The black lines show the 
sum of all three components, and the dotted lines indicate 
the ratio required for perfect virialization. In the minihalo 
case, the total energy is dominated by thermal energy, al- 
though its share decreases towards the centre where cooling 
via molecular hydrogen becomes important. The radial ki- 
netic energy dominates over the tangential component down 
to r tan ~ 5 pc, where the mass fraction of cold gas rapidly 
rises and the cloud becomes rotationally supported. Effi- 
cient cooling implies that the total energy drops below that 
required for perfect virialization. 

In the atomic cooling halo, the total energy at r v - u — 
1 kpc is dominated by bulk radial motions. The distinction 
between hot and cold gas in the right panels of Fig. 11 shows 
that a large fraction of the kinetic energy injected into the 
galaxy comes from cold gas accreted along filaments, even 
though its mass fraction is initially small. The energy in tan- 
gential motions begins to dominate at r tan — 200 pc, show- 
ing that the radial energy of the cold gas flowing in along 
filaments is converted into turbulent motions. This is fun- 
damentally different from the collapse of gas in minihaloes, 
where the radial energy is converted into a directed rota- 
tion along a single axis. The distinct features at r ~ 350 pc 
are caused by a subhalo that has not yet merged with the 
central clump (see also Figs 2 and 3). Once again, the total 
energy budget falls below that required for perfect virializa- 
tion, as atomic hydrogen as well as molecular cooling are 
able to radiate away a significant fraction of the potential 
energy released. We conclude that the high energy input by 
cold accretion is ideally suited to drive turbulence at the 
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Figure 8. The central ~ 40 kpc (comoving) of the computational 
box, roughly delineated by the insets in Fig. 2. Shown is the 
Mach number in a slice centred on the BH at the centre of the 
galaxy, indicated by the filled black circle. The dashed line denotes 
the virial radius at a distance of ~ 1 kpc. The Mach number 
approaches unity at the virial shock, where gas accreted from the 
IGM is heated to the virial temperature over a comparatively 
small distance. Inflows of cold gas along filaments are supersonic 
by a factor of ~ 10 and generate a high amount of turbulence at 
the centre of the galaxy, where typical Mach numbers are between 
1 and 5. 

centre of the galaxy, where bulk radial inflows are converted 
into turbulent motions on small scales. 

5.2 Shocks and fragmentation properties 

Shock fronts can arise where supersonic flows experience 
sudden deceleration and, if unorganized, indicate the pres- 
ence of supersonic turbulence. As discussed above, cold ac- 
cretion is a viable agent for driving turbulence, due to 
the prodigious amount of momentum and kinetic energy it 
brings to the centre of the galaxy. In Fig. 12, we show the 
divergence and vorticity of the velocity field during the viri- 
alization of the galaxy. A comparison with Fig. 8 implies 
that there are indeed regions of supersonic flow that expe- 
rience rapid deceleration and form shocks. In our case, two 
physically distinct mechanisms are responsible for creating 
these shocks. The virial shock forms where the ratio of in- 
fall velocity to local sound speed approaches unity, and is 
clearly visible in the left-hand panel of Fig. 12. The veloc- 
ity divergence is negative since the gas rapidly decelerates, 
while the vorticity is almost negligible. In contrast, the unor- 
ganized multitude of shocks that form near the centre of the 
galaxy are mostly caused by accretion of cold, high-velocity 
gas from filaments. These are more pronounced than the 
virial shock and have a significantly higher angular com- 
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S\zt: 40 Mpc {removing) 
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z - 10,62 
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Figure 9. The central ~ 40 kpc (comoving) of the computational box, roughly delineated by the insets in Fig. 2. Shown is the hydrogen 
number density (left-hand panel) and temperature (right-hand panel) in a slice centred on the BH at the centre of the galaxy, indicated 
by the filled black circle. The dashed lines denote the virial radius at a distance of ~ 1 kpc. Hot accretion dominates where gas is accreted 
directly from the IGM and shock-heated to ~ 10 4 K. In contrast, cold accretion becomes important as soon as gas cools in filaments and 
flows towards the centre of the galaxy, such as the streams coming from the left- and right-hand side. They drive a prodigious amount 
of turbulence and create transitory density perturbations that could in principle become Jeans-unstable. In contrast to minihaloes, 
the initial conditions for second-generation star formation are highly complex, with turbulent velocity fields setting the fragmentation 
properties of the gas. 



ponent. They create transitory density perturbations that 
could in principle become Jeans-unstable and trigger the 
gravitational collapse of individual clumps. 

How does the turbulence generated in the infalling ma- 
terial influence its fragmentation behaviour and control sub- 
sequent star formation? From detailed observational and 
theoretical studies of star formation in our Milky Way we 
know that turbulence plays a pivotal role in the formation 
of stars and star clusters. It is usually strong enough to 
counterbalance gravity on global scales. By the same token, 
however, it will usually provoke collapse locally. Turbulence 
establishes a complex network of interacting shocks, where 
regions of high density build up at the stagnation points 
of convergent flows. To result in the formation of stars, lo- 
cal collapse must progress to high enough densities on time 
scales shorter than the typical interval between two succes- 
sive shock passages. Only then can the collapsing core de- 
couple from the ambient flow pattern and build up a star. 
The accretion flow on to these objects and consequently the 
final stellar mass strongly depends on the properties of the 
surrounding turbulent flow. 

In concert with the thermodynamic properties of the 
gas, leading to the cooling of high-density material to the 
CMB limit (see Section 4), length scale and strength of the 
turbulence are the most important parameters governing 
its fragmentation behaviour and consequently the proper- 
ties of star formation, such as its timescale and overall ef- 



ficiency (Klesse n et al, Il200d : IVazquez-Semadeni et af1l2003l ; 
iKrumholz fc McKedbOOoT T In the atomic cooling halo dis- 
cussed here, this will eventually lead to the transition to 
Pop II star formation. However, a quantitative understand- 
ing of the fragmentation behaviour of the turbulent gas 
would require dedicated high-resolution simulations, which 
is beyond the scope of this work. 



6 MASSIVE BLACK HOLE GROWTH 

Galaxy formation in general involves the co-evolution of a 
central black hole and the surrounding stellar system, the 
one influencing the other. Two crucial unsolved problems 
are: What were the seeds for BH growth, and how impor- 
tant was this co-evolution at very high redshifts? We here 
begin to address these questions. Different scenarios have 
been suggested to account for the seeds of BH growth (|Reesl 
|l984| ): the direct collapse of gas in atomic cooling haloes 
in the presence of a strong photodissociating background, 
or stellar remnants of massive, metal-free stars. In the fol- 
lowing, we discuss the latter possibility and investigate the 
growth of a MBH forming at the centre of the galaxy. 

6.1 Accretion rate 

Even though studies of stellar evolution have shown that 
primordial stars may explode as energetic SNe, we here 
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Figure 10. The central ~ 40 kpc (comoving) of the computational box, roughly delineated by the insets in Fig. 2. Shown is the radial 
(left-hand panel) and tangential velocity in the x-y plane (right-hand panel) in a slice centred on the BH at the centre of the galaxy, 
indicated by the filled black circle. The dashed lines denote the virial radius at a distance of ~ 1 kpc. Streams of cold gas from filaments, 
such as those coming from the left- and right-hand side, are clearly visible and can have velocities of up to 20 km s — x . Some even 
penetrate the centre and create regions of positive radial velocities. Angular velocities are particularly high towards the centre of the 
galaxy, where bulk radial inflows are converted into turbulent motions on small scales. The presence of flows in both directions implies 
that these are unorganized instead of coherently rotating, such as is the case in minihaloes (see also Fig. 11). 



assume that all Pop H I T stars collapse d irectly to BHs 
jHeger fc Wooslevl |2002| ; iHeger et all hoosl ). Their initial 
mass is dictated by the resolution limit to Mbh — 2 x 
10 3 Mq , with accretion on to the BH governed by the crite- 
ria discussed in Section 2. Recent investigations have shown 
that photoheating by the progenitor star can delay effi- 
cient accretion by reducin g the central density to < 1 cm -3 
I Johnson fc Bromnj 12007m . However, the suppression of ac- 
cretion also depends on the detailed merger history of the 
host halo. For example, a major merger occuring just after 
the formation of the BH could transport enough cold gas 
to its centre to enable accretion at the Eddington rate. As 
shown in Fig. 6, some minihaloes merge shortly after form- 
ing a BH, in some cases after only a few million years. On 
the other hand, if mergers are absent or incoming haloes are 
not sufficien tly massive, accretion cou ld be suppressed for 
> 100 Myr l| Johnson fc Bromml |2007m . In our approach to 
derive an upper limit on the accretion rate, we neglect the 
effects of photoheating by the progenitor star such that ac- 
cretion is governed solely by the supply of cold gas brought 
to the centre of the halo. 

A more precise modelling would also require a prescrip- 
tion for radiation emitted by the BH-po wered miniquasar 
and i ts feedback on the surrounding disc (|Volonteri fe Reesl 
2006). In a simpler approach, one can assume a given radia- 
tive efficiency e, which denotes the ratio of BH luminosity to 
accreted mass energy, and assume that accretion is spheri- 
cally symmetric. This leads to the Eddington accretion rate: 



• _ 1 M BH 

AiEdd = — , 

£ ISalp 

where tsaip is the Salpeter time, defined by: 

CO"T 



^Salp = 



450 Myr 



(3) 



(4) 
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The mass of the BH as a function of time is thus given by 



Mbh (t) = M BH (to) exp 



1 - e t - t 

e tsalp 



(5) 



In Fig. 13, we compare the mass growth of the most mas- 
sive BH with the Eddington-limited model, using a fiducial 
value of e = 1/10. We find that the accretion rate remains 
roughly constant at ~ 5 x 10~ 3 M Q yr" 1 , such that the 
BH grows from ~ 2 x 10 3 to ~ 10 6 M in the course of 
~ 300 Myr. Due to our neglect of radiative feedback, the BH 
accretes well above the Eddington rate throughout most of 
its lifetime. At later times the accretion rate stagnates, most 
likely caused by the high kinetic energy input at the centre 
via cold accretion (see Fig. 11). Consequently, the fraction 
of unbound gas near the BH increases and its mass growth 
is slowed. We conclude that Eddington accretion is possi- 
ble under the most favourable circumstances, for example, 
where a recent merger brings an ample supply of cold gas 
to the centre of the halo, but generally radiative feedback 
by the progenitor star and the dis c around the BH will lead 
to sub-Eddington acc retion rates l|johnson fc Brommll2007l : 
iPelupessv et al.ll2007l ). 
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Figure 11. The energy distribution and mass fraction of cold (< 500 K) versus hot (> 500 K) gas in radial shells for the first star-forming 
minihalo just before the formation of the sink particle, and the atomic cooling halo assembling at z ~ 10. The azimuthally averaged 
ratio of radial, tangential and thermal to potential energy are shown as blue, green and red lines, respectively. The black lines show the 
sum of all three components, and the dotted lines indicate the ratio required for perfect virialization. In the minihalo case, the total 
energy is dominated by thermal energy, although its share decreases towards the centre where cooling via molecular hydrogen becomes 
important. The radial kinetic energy dominates over the tangential component down to rt an — 5 pc, where the mass fraction of cold gas 
rapidly rises and the cloud becomes rotationally supported. In the atomic cooling halo, the total energy at r v i r ~ 1 kpc is dominated 
by bulk radial motions. The distinction between hot and cold gas in the right panels shows that most of the kinetic energy injected into 
the galaxy comes from the cold gas accreted along filaments, even though its mass fraction is initially small. The tangential component 
begins to dominate at rt a n — 200 pc, where the radial flow of cold gas along filaments is converted into turbulent motions. The distinct 
features at r ~ 350 pc are caused by a subhalo that has not yet merged with the central clump (see also Figs 2 and 3). 



6.2 Accretion luminosity 

The radiation generated by the BH-powered miniquasar can 
have numerous effects on the formation of the first galax- 
ies. In addition to its negative feedback on star forma- 
tion via photoheating, the emitted radiation can contribute 
to the LW b ackground, as well as to the reionization of 
the Universe llMadau fc Quataertj 120041 ; iRicotti fc Ostrikerl 
12004 iKuhlen fc Madaull2005t ). In the following, we derive 
an upper limit on the photodissociating flux and the num- 
ber of ionizing photons emitted by a stellar remnant BH 
accreting at the Eddington limit. 



6.2.1 Photodissociating flux 

To determine the flux of LW photons, we first model the 
temperature profile of the surrounding accretion disc: 



T(r) = 



3 GMbhMbh 
8tt 



(6) 



where r is the distance from the BH, Mbh its mass, Mbh 
the accretion rate and ctsb the Stefan-Boltzmann constant 
|Pringlelll98"ll ). For simplicity, we have taken the disc to be 



a thin disc, such that each annulus radiates as a blackbody 
of temperature given by the above equation. The inner-most 
radius of the disc is given by: 



2 km 



-Me 



M . 



(7) 



correspon ding to a high value for the BH spin paramete r 
a > 0.9 (|Makishima et all |2000| : IVierdavanti et all I2007T I. 
which we expect considering the large angular momentum 
of the accreted gas (see Fig. 11). We integrate the flux over 
the surface of the disc from TWer to r out cr = 10 4 rj nne r, 
where the contributions to both the photodissociating and 
ionizing fluxes are negligible. To determine the flux emit- 
ted in the LW bands, we evaluate the total emitted flux 
at 12.87 eV. In the upper panel of Fig. 14, we show the 
LW flux in units of 10~ 21 erg s _1 cm -2 Hz -1 sr _1 for 
the case of Eddington-limited accretion at 1 kpc distance 
from the BH. We consider initial BH masses of 100 and 
500 Mq, roughly the range expe cted for the direct col- 
lapse of massive Pop III stars l|Heger fc Wooslevl 120021 : 
iHeger et all 120031 ). As Fig. 14 shows, J LW can greatly ex- 
ceed the critical value of > 10~ 2 required for the suppres- 
sion of star formation in minihaloes, which mostly relies 
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Figure 12. The central ~ 40 kpc (comoving) of the computational box, roughly delineated by the insets in Fig. 2. We show the divergence 
(left-hand panel) and z-component of the vorticity (right-hand panel) in a slice centred on the BH at the centre of the galaxy, indicated 
by the filled black circle. The dashed lines denote the virial radius at a distance of ~ 1 kpc. The most pronounced feature in the left-hand 
panel is the virial shock, where the ratio of infall speed to local sound speed approaches unity and the gas decelerates over a comparatively 
small distance. In contrast, the vorticity at the virial shock is almost negligible. The high velocity gradients at the centre of the galaxy 
indicate the formation of a multitude of shocks where the bulk radial flows of filaments are converted into turbulent motions on small 
scales. 



on efficient H2 cooling IpVIachacek et al. 200 ll; Yoshida et al.l 
120031 : lO'Shea fc Normanll200Sl : I Johnson et alll2007bh . Fur- 
thermore, even a modest LW flux can dissociate enough H2 
such that HD formation and cooling never become impor- 
tant, reducing the t emperature to which the gas can cool 
I Yoshida et all l2007h . The impact on primordial star for- 
mation in BH-hosting galaxies might thus be severe even 
for sub-Eddington accretion rates, implying that this effect 
must be taken into account in future work. 



Another important issue concerns the contribution to 
the global LW background. We may estimate a maximum 
global LW background by assuming that each atomic cool- 
ing halo at 2 > 10 hosts a BH accreting at the Eddington 
rate. We can then find an upper limit to the global LW back- 
ground by summing up the contributions from BHs within 
a distance equal to the maximum mean free path for a LW 
p hoton, r max ~ 1 Mpc a t z = 15, following the prescription 
in Ijohnson et all l|200 7bl . This estimate yields a maximum 
global LW background comparable to the LW flux from a 
single source, as shown in Fig. 14. Such a high flux could 
have profound consequences for further star formation in 
minihaloes, but in most cases radiative feedback by the pro- 
genitor star will significantly delay accretion, such that a 
global LW background fuelled by accretion on to BHs will 
likely be subdominant compared to a stellar LW background 
l|Pelupessv et alj|2007h . 



6.2.2 Ionizing flux 

The amount of ionizing radiation released by the accret- 
ing BH can be determined in analogy to our calculation of 
the LW flux. An integration over the temperature profile 
of the accretion disc yields the total number of hydrogen- 
ionizing photons emitted per second, shown in the bottom 
panel of Fig. 14. While massive Pop HI stars emit of the 
order of 10 50 hydrogen-ionizi ng photons per second, these 
stars live for only < 3 Myr (Br omm et all 1200 ll ; ISchaereil 
2002). However, as Fig. 14 shows, if Pop III relic BHs are 
able to accrete efficiently, they may emit 10 - 100 times 
this number for > 100 Myr. This enormous flux of ion- 
izing radiation could power H II regions with radii of the 
order of 10 kpc, larger and longer-live d than the transient 
H ll regions of individ ual Pop III stars (Yo shida et alll2007l ; 
I Johnson et~alll2007af ). Star formation in minihaloes within 
the H ll region could be suppressed if accretion is continuous 
and drives a persist ent radiative flux (|Ahn fc Shapirdl2007l ; 
IWhalen et alll2007h . Due to this dramatic radiative feedback 
associated with high accretion rates, though, the gas in the 
protogalaxy is likely to be heated and driven away from the 
BH, once again result ing in sub-Eddington accretion rates 
l|Pelupessv et alll2007h . 

We have also calculated the number of He n-ionizing 
photons emitted by the accreting BH, and find that this is 
within a factor of < 2 of the number of hydrogen-ionizing 
photons, owing to the high temperatures of the accretion 
disc. Thus, if Pop III relic BHs are able to accrete efficiently, 
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Figure 13. The mass (top panel) and accretion rate (bottom 
panel) of the central BH as a function of redshift, shown for the 
simulation (solid lines) and Eddington-limited accretion (dotted 
lines). The accretion rate remains roughly constant at ~ 5 X 
10~ 3 Mq yr _1 , such that the BH grows from ~ 2 X 10 3 to ~ 
10 6 Mq in the course of ~ 300 Myr. This is a strict upper limit as 
radiation effects are not taken into account. Accretion is initially 
super-Eddington due to the high amount of cold gas brought to 
the centre of the galaxy, but stagnates once turbulence becomes 
important and the fraction of unbound gas increases. 

they may also contribute to the reionization of helium, driv- 
ing He hi regions that can be as large as their H II regions 
ijFurlanetto fc Ohll2007f l. 



7 SUMMARY AND CONCLUSIONS 

We have investigated the formation of the first galaxies 
with highly resolved numerical simulations, taking into ac- 
count all relevant primordial chemistry and cooling. The first 
galaxies form at redshifts z > 10 and are characterized by 
the onset of atomic hydrogen cooling, once the virial tem- 
perature exceeds ~ 10 4 K, and their ability to retain pho- 
toheated gas. We have described the merger history of a 
~ 5 x 10 r Mq system in great detail and found that in the 
absence of stellar feedback 10 Pop III. 1 stars form in mini- 
haloes prior to the assembly of the galaxy. Infalling gas is 
partially ionized at the virial shock and forms a high amount 
of H2 and HD, allowing the gas to cool to the temperature 
of the CMB and likely form Pop III. 2 stars with > 10 Mq. 
Accretion on to the galaxy proceeds initially via hot ac- 
cretion, where gas is accreted directly from the IGM and 
shock-heated to the virial temperature, but is quickly ac- 
companied by a phase of cold accretion, where the gas cools 
in filaments before flowing into the parent halo with high ve- 
locities. The latter drives supersonic turbulence at the centre 
of the galaxy and thus plays a key role for second-generation 
star formation. Finally, we have investigated the growth of 
BHs seeded by the stellar remnants of Pop III. 1 stars and 
found that accretion at the Eddington limit might be pos- 
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Figure 14. The radiation due to thermal emission from the ac- 
cretion disc surrounding a Pop III relic BH accreting at the Ed- 
dington limit, for an inital mass of 100 (solid lines) and 500 Mq 
(dotted lines). Top panel: The photodissociating flux, Jlw 5 in 
units of 10 -21 erg s — 1 cm -2 Hz -1 sr" 1 , at a distance of 1 kpc 
from the BH. Bottom panel: The number of hydrogen-ionizing 
photons emitted per second from the accretion disc, in units of 
10 50 s- 1 . 

sible under the most favourable circumstances, but in most 
cases radiation emitted by the progenitor star and the accre- 
tion disc around the BH will lead to sub-Eddington accretion 
rates. 

Depending on the strength of radiative and SN-driven 
feedback, some galaxies might remain metal-free and form 
intermediate-mass Pop III. 2 stars. The inclusion of radia- 
tive feedback would likely increase the fraction of Pop III. 2 
material, as it enhances the degree of ioniza tion, and, con- 
sequen t ly, the amount of m olecules formed l|johnson et al.l 
l2007al ; IWise fc Abel l2007d ). Observational signatures of 
intermediate-mass primordial stars might include gamma- 
ray bursts ( GRBs) for the case of a rapidly rotating pro- 
genitor (e.g. iBromm fc Loebl 2006). or distinct abundance 
patte rns produced by core-co llapse SNe experiencing fall- 
back l|Umeda fc NomotollioO^ ). However, since the number 
of Pop III. 1 stars formed in minihaloes prior to the assem- 
bly of the galaxy is of the order of 10, it seems much more 
likely that at least one star will end in a violen t SN explosion 
and pre-enrich the halo to supercritical leve ls (|Bromm et al.l 
120051 ; iGreif et al.ll2007l ; FWise fc Abelll2007cl ). In combination 
with the onset of turbulence, metal mixing in the first galax- 
ies will likely be highly efficient and co uld lead to the for - 
mation of the first low-mass star clusters l|Clark et al.ll2008h , 
in ex treme cases possibly eve n to metal-poor globular clus- 
ters (|Bromm fc Clarke] I2002T I. Some of the extremely iron- 
deficient, but carbon and oxygen-enhanced stars observed in 
the halo of the Milky Way may thus have formed as early 
as redshift z ~ 10. 

In future work, we plan to include the effects of ra- 
diative and SN-driven feedback by previous star formation 
in minihaloes, as well as the distribution of metals and its 



© 0000 RAS, MNRAS 000, 000-000 



The First Galaxies 17 



concomitant cooling. We will study the fragmentation of 
gas accumulating at the centre of the galaxy and address 
the issue of turbulence-driven star formation in detail. The 
goal of making realistic predictions for the first generation 
of starburst-galaxies, to be observed with the James Webb 
Space Telescope (JWST), is clearly coming within reach. 
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